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Abstract Antiproton annihilations on nuclei provide a very interesting way to 
study the behaviour of strange particles in the nuclear medium. In low energy p 
annihilations, the hyperons are produced mostly by strangeness exchange mecha- 
nisms. Thus, hyperon production in pA interactions is very sensitive to the prop- 
erties of the antikaon-nucleon interaction in nuclear medium. Within the Giessen 
Boltzmann-Uehling-Uhlenbeck transport model (GiBUU), we analyse the experi- 
mental data on A and K s production in pA collisions at pi a b =0.2 — 4 GeV/c. A 
satisfactory overall agreement is reached, except for the Kg production in p+ 20 Ne 
collisions at pi a b = 608 MeV/c, where we obtain substantially larger Kg pro- 
duction rate. We also study the 3 hyperon production, important in view of the 
forthcoming experiments at FAIR and J-PARC. 

Keywords pA collisions • A, Kg and 3 production • BUU model 



1 Introduction 

About two decades ago several experiments have been done on strangeness produc- 
tion in p-nucleus reactions. In the hydrogen bubble chamber experiment at BNL 
[T], the yl-hyperon production in collisions p(0 - 450MeV/c) + 12 C, 48 Ti, 181 Ta and 
208 Pb has been measured. At KEK 2,3|, the hydrogen bubble chamber measure- 
ments have been performed for the Kg, A and A production from p(4 GeV/c) 
+ 181 Ta. Finally, at LEAR _4 , the K% and A rapidity yields from p(607 MeV/c) + 
20 Ne interactions have been studied using the streamer chamber filled with a nat- 
ural 20 Ne gas. Recently, K production in annihilation of stopped p's on p, d, 
3 He and He has been measured at LEAR [5 applying the magnetic spectrometer 
Obelix. 
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These experiments revealed some interesting features which still remain to be 
explained by theory: 

— Large ratio of the particle yields A/K s = 2 — 3 both for light ( 20 Ne) and heavy 
( 181 Ta) targets. 

— A rapidity spectrum peaked close to the target rapidity even for energetic 
collisions p(4 GeV/c) + 181 Ta. 

— Enhanced strangeness production for p annihilations at rest involving more 
than one nucleon (baryon number of annihilating system B > 0) . 

In order to describe the yields of strange particles and, in particular, the large 
A/Kg ratio in p + 181 Ta interactions at 4 GeV/c, Rafelski [6] assumed an anni- 
hilation fireball in the state of a supercooled quark-gluon plasma (QGP) propa- 
gating through the target nucleus and absorbing nucleons. Eventually, this fireball 
reaches a high baryon number, B ~ 10, and then hadronizes producing the excess 
of strange particles. Such an exotic scenario has been questioned in the following- 
up theoretical work by Cugnon et al. [7], where the intranuclear cascade (INC) 
calculations of the strangeness production in p-nucleus interactions have been per- 
formed. 

In this talk, we present the results of our transport-theoretical analysis of ex- 
perimental data on strangeness production from p-nucleus interactions in-flight. 
We also make predictions on the S-hyperon (S = —2) and yl-hypernuclear pro- 
duction. Section [2] contains a brief description of the GiBUU model applied in our 
calculations. In Section [3] we present the numerical results and discuss them. The 
summary is given in Section [4] 

2 Model 

The GiBUU transport model [8 solves the coupled system of relativistic kinetic 
equations for the different sorts of hadrons: 

(p5) _1 [pX + (p*^T + m*(d>*))£n/i(z, P *) = /,[{/}] , (i) 

where fi(x, p*) with x = (i, r) denotes the distribution function of the particles 
of sort i (i = N, N, A, A, Y, Y, S, S, it, 77, lu, p, K, K,... taking into 
account isospin projections) in the six-dimensional phase space (r, p*). If the r.h.s. 
of Eq. (fl|), i.e. the collision term U [{/}], would be zero, then Eq. (fTJ) would be a 
classical Vlasov equation describing the propagation of the particles in the mean 
field potentials of a nuclear and electromagnetic nature. 

The relativistic form of Eq. ([TJ) needs to explain some specific quantities: p* M = 
p M — Vf 1 is the kinetic four-momentum, F? v = d^V" — d v Vf i is the field tensor, 
and m* = mi + Si is the effective mass. The particles are assumed to be on their 
in-medium mass shells, p^p*^ = m* 2 . The scalar and vector fields are expressed, 
respectively, as Si = g a io and = gui^ + gpi^p 3 ^ + qiA^. Here, a (I = 0, J = 
0), w M (I = 0, J = 1) and (I = 1, J = 1) are the mean mesonic fields and 
A^ = (A°,0) is the Coulomb field. 

The mean mesonic and Coulomb fields are calculated, respectively, from the 
static Klein-Gordon-like and Poisson equations with the source terms provided 
by the space- and time-dependent particle densities and currents. The meson- 
nucleon coupling constants and the self-interaction parameters of the cr-field are 
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adopted from the NL3 version of the non-linear Walecka model [9,. The coupling 
constants of the other baryons with mean mesonic fields are obtained from the 
corresponding meson-nucleon coupling constants by simple rescaling taking into 
account the light-quark contents, and, for antibaryons, the G-parity symmetry 
and phenomenological depth of the antiproton optical potential (Re(V r op t) — —150 
MeV in the nuclear center). This leads, e.g., to a A(A) potential of ~ — 40(— 450) 
MeV and K(K) potential of ~ -20(-220) MeV (see refs. [T0II11II12] for further 
details of the mean fields) . 

The collision integral in the r.h.s. of Eq. (fl|) describes the antibaryon-baryon an- 
nihilation, elastic and inelastic hadron-hadron scattering processes, and resonance 
decays. It includes the elementary hadron-hadron cross sections and resonance 
widths as parameters. The elementary cross sections are obtained by phenomeno- 
logical fits to experimental data or given by theoretical calculations. 

The following antibaryon-baryon reaction channels are included: BB — » mesons 
simulated by statistical annihilation model [13] . BB — > BB (elastic and charge 
exchange), NN <->• NA(AN), NN -> A A, NN{NA, AN) -> AE(SA), NN(NA, 
AN) SB. For invariant energies y/s > 2.4 GeV (pi a b > 1-9 GeV/c for NN), 
the inelastic production in antibaryon-baryon collisions, B\Bi — > B^Bi + mesons, 
is simulated with the help of the FRITIOF model. 

The baryon-baryon reactions implemented in the model are: BB — » BB (elastic 
and charge exchange), NN <H> NNir, NN o AA, NN o NR, where R denotes 
any nonstrange baryonic resonance, N(A,N*)N(A,N*) -> N(A)YK, YN -> 
YN, BN -)■ AA, BN ->■ AS, BN ->■ SN. For > 2.4 GeV, the inelastic 
production in baryon-baryon collisions, B1B2 — > B3B4 + mesons, is simulated by 
using the PYTHIA model. 

The meson-baryon collisions taken into account in the model are: ttN o R, 
nN KKN, Tr(n,p,uj)N YK, KN O Y* , KN -> KN, KN O Ytt, KN O 
Y*ir, KN —> BK. For ^fs > 2.2 GeV, the meson-baryon collisions are simulated 
by applying the PYTHIA model. 

The GiBUU model also includes the strangeness production/absorption chan- 
nels in meson-meson collisions: MM <-> K(K*)K(K*), where M denotes any 
nonstrange meson (M = 7r, r\, 77', a, p, uj,). Further details on the elementary 
cross sections included in the model can be found on the GiBUU web-site [8] and 
in the review article [ll] . 

3 Results 

We start with the low beam momenta and consider the annihilations of p at 225 
MeV/c on several target nuclei. Figure [1] shows the probability of the A production 
per p annihilation event on a nucleus as a function of the nucleus mass number. 
We have separated the A hyperons emitted to free space from those bound in 
the nuclear remnant at the end of time evolution (~ 200 fm/c starting from p 
at the longitudinal distance of nuclear radius + 5 fm from the nuclear centre). 
Our results on free yl's agree quite well with experiment. The probability of the 
bound A- (^1^1-) nuclear system production, as expected, grows with the target 
mass number reaching ~ 3% (0.01%) for the lead target. 

In Fig.[2l we present the rapidity spectra of (A+S°) and Kg for p annihilations 
on 20 Ne at 607 MeV/c. Since the S° — > A-y decay is not taken into account 
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p(225 MeV/c)+C, Ti, Ta, Pb -> A(2A) X 
10 -1 , ^ 
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Fig. 1 The cross sections of the free A production, bound A production, and two bound A'b 
production for p(225 MeV/c)+ 12 C, 48 Ti, 181 Ta and 208 Pb. The cross sections are normalized 
on the annihilation cross sections of p on the corresponding nuclei. The data for free A's pLj 
are for the beam momentum range 0-450 MeV/c. 



in GiBUU, we have summed up both A and S spectra in order to compare 
with experimental data. The data on /1-rapidity yields are described quite well by 
GiBUU. By studying the Y(Y*) production rate more in-detail, we have found 
that ~ 80% of this rate is due to the strangeness exchange reactions of the type 
KN ->• Yn. 

The calculated K$ yield, N K o = (N K o + Ng- )/2, is considerably higher than 
experiment, as one can see from Fig. [2] The INC calculations [7] are also in a rather 
good agreement with A rapidity spectrum, but overpredict the Kg production. 

Fig. E] presents the (A + K% and (A + S°) rapidity spectra from p+ 181 Ta 
collisions at 4 GeV/c. The GiBUU model reproduces the tails of the y-distribution 
of A's fairly well, but underestimates the data at the peak position, y ~ 0.3. On 
the other hand, the INC model overestimates the data around the target rapidity, 
y = 0. Both models produce the peak position of the A rapidity spectrum at 
y ~ 0. The theoretical rapidity spectra of Kg agree with data reasonably well, 
except for target rapidities, where some excess of Kg is still visible. By inspecting 
the Y(Y*) production rate we again found that ~ 70 — 80% of this rate is due to 
the strangeness exchange processes K(K*)B —> YX, KB — > Y* , and KB — > Y*n, 
while the total contribution of the BB collision channels — including the direct 
NN channel — is on the level of a few percent only. The difference between the 
(A+S ) rapidity spectra produced by GiBUU and INC is mostly due to somewhat 
different angular distributions for the NN — >• YY processes and due to different 
YN annihilation cross sections used in the both models. We also note that the 
INC model does not contain mean field potentials, while our calculations include 
strongly attractive antihyperon potentials. 
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Fig. 2 The rapidity distributions of the (A+H°) hyperons and K% from p(607 MeV/c)+ 20 Ne 
annihilations. The distributions are normalized so that their integrals over y give the number 
of particles per annihilation event. The INC calculations from [7] are also shown. The data 
(histograms) are from ref. [4]. Vertical arrows indicate the rapidity of a p-nucleon center-of- 
mass (cm.) system. 



Let us, finally, discuss the predictions of the GiBUU model for the S = —2 
hyperon production in p+ 197 Au collisions at 3 GeV/c. We have found that the pre- 
dominant channel of H-hyperon production is the collisions of the strange mesons 
K, K* , K and K* with baryons which are responsible for ~ 35% of the to- 
tal S yield. Almost of the same importance are the S* —> decays (~ 26%), 
and strange meson-hyperon collisions (~ 17%). The direct channel NN —> SS is 
responsible only for a quite small fraction of produced S'a (~ 5%). 
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Fig. 3 The (A + E°), K° s and (A + E°) rapidity distributions for the collisions 
p(4 GeV/c) 181 Ta. The INC calculations are taken from [7]. The data are from [3]. 
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Fig. 4 Rapidity distributions of (A + 17°), £ and from p interactions with 197 Au at 3 
GeV/c. The spectrum of S~ hyperons is multiplied by a factor of 10. 



Fig-H]shows the rapidity spectra of S~ , (A+S°) and Kg for the p(3 GeV/c) + 
197 Au collisions. The spectra of the (A + S°) hyperons are peaked close to y = 0. 
The Kg spectra have also a maximum at y ~ spread towards forward rapidities. 
This behaviour has been interpreted in [6] as a signature of strangeness production 
from a QGP hreball. In this picture, the rapidity spectra of all strange particles 
emitted from the fireball should have their maxima at the cm. rapidity of the 
fireball. However, the spectrum of S~ hyperons produced by GiBUU calculations 
is shifted to the forward rapidities and peaked at y ~ 0.6. The difference between 
the (A + S°) and E~ rapidity spectra can be understood in our pure hadronic 
transport calculations as the consequence of different thresholds for Y and a 
production in antikaon-nucleon collisions. Indeed, the S = — 1 hyperons are mostly 
produced in exothermic reactions, like KN — > Y-n, with slow initial K. As a result, 
the produced hyperon Y is also slow and its momentum is isotropically distributed 
in the laboratory frame. However, the double strangeness exchange process, KN — > 
SK, mainly responsible for S production, is endothermic (pf^b = 1-05 GeV/c — 
K beam momentum at threshold, 2/c.m. = 0.55 — cm. rapidity at threshold) and 
requires a fast incoming K. Such antikaons are mostly emitted in the forward 
direction in the laboratory frame. Thus, the outgoing S~ moves also forward due 
to the cm. motion of the KN system. 
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4 Summary 

We have performed microscopic transport calculations of strangeness production 
in antiproton-nucleus collisions at pi a b = 0.2 — 4 GeV/c on the basis of the GiBUU 
model. The main results can be summarized as follows: 

— The experimental data on /1-yields are reasonably well described. However, 
there are some deviations in detailed shape of yl-rapidity spectrum at 4 GeV/c. 
The i£g-yields are overestimated. The yield ratios for p+ 20 Ne at 608 MeV/c 
are: A/K% = 1.0 (GiBUU), A/K% = 1.4 (INC), A/K° s = 2.3 ±0.7 (exp.). This 
indicates missing antikaon absorption by KN — » Yn in transport calculations. 

— The peak positions of the (A + S )- and H-hyperon rapidity spectra strongly 
differ. In future experiments at FAIR and J-PARC, this can be tested in order 
to support or exclude the possible exotic mechanism of strangeness production 
via a QGP fireball. 

— A- and /1/1-hypernuclei production is possible on the primary target at low 
p-beam momenta. 
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